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ABSTRACT 

14 The flybys of Jupiter by the Voyager spacecraft in 1979, and over two decades later 

15 by Cassini in 2000, have provided us with unique datasets from two different epochs, 

16 allowing the investigation of seasonal change in the atmosphere. In this paper we model 

17 zonal averages of thermal infrared spectra from the two instruments. Voyager 1 IRIS 

18 and Cassini CIRS, to retrieve the vertical and meridional profiles of temperature, and 

19 the abundances of the two minor hydrocarbons, acetylene (C2H2) and ethane (C2H6). 

20 The spatial variation of these gases is controlled by both chemistry and dynamics, and 

21 therefore their observed distribution gives us an insight into both processes. We find that 

22 the two gases paint quite different pictures of seasonal change. Whilst the 2-D cross- 

23 section of C2H6 abundance is slightly increased and more symmetric in 2000 (northern 

24 summer solstice) compared to 1979 (northern fall equinox), the major trend of equator 

25 to pole increase remains. For C2H2 on the other hand, the Voyager epoch exhibits almost 

26 no latitudinal variation, whilst the Cassini era shows a marked decrease polewards in 

27 both hemispheres. At the present time, these experimental findings are in advance of 

28 interpretation, as there are no published models of 2-D Jovian seasonal chemical variation 

29 available for comparison. 

key words=JUPITER ATMOSPHERE; ATMOSPHERIC ABUNDANCES, 
OUTER PLANETS; INFRARED SPECTROSCOPY; ABUNDANCE RETRIEVAL 



30 1 Introduction 

31 The Voyager 1 and 2 encounters with Jupiter of March and July 1979 (jovian northern 

32 fall equinox) offered the first opportunities to map Jupiter in the thermal infrared with 

33 a latitude resolution of better than 10° - then unobtainable from the ground - and a 

34 spectral resolution sufficient to measure the abundances of trace gas spec ies. The too 

35 f or th is mapping was IRIS (the InfraRed Interferometer and Spectrometer, iHanel et al. 



19771 ). and trace gases of interest included ethane (C2H6) and acetylene (C2H2) - secondary 



37 species derived from photolysis of the primary carbon-bearing molecule, methane (CHd) . 



38 Altho ugh some preliminary findings were published in a conference report (JMaguire et al. 



19841 ). a full radiative transfer model was never applied to retrieve a meridional trend in 



40 the abundances. 

41 Two decades later, the Cassini spacecraft swung by Jupiter in a distant fiyby maneuver 

42 (137 Rj) en route to Saturn. The official encounter period, during which time observations 

43 were made, lasted six months symmetric about the closest approach on December 30th 

44 2000 (northern summer solstice). Carrying on board the Composite Infrared Spectrometer 

45 (CIRS), Cassini was able to re- map the planet in the thermal infrared, achieving a spatial 

46 resolution comparable to Voyager IRIS but with much higher maximum spectral resolution 

47 (0.48 cm~^, versus 3.9 cm~^ for IRIS, full width to half maximum (FWHM)). Meridional 

48 trends of stratospheric and tropospheric C2H2 and C2H6 have been r etrieved from CIRS 



Nixon et al. 



(120071 ) 



49 data acquired during December 1-15 2000 prior to cloest approach, by 

50 (hereafter Paper I). 

51 The abundances of the minor hydrocarbons are important because these species serve 
62 as tracers of atmospheric circulation. The photochemical lifetime of ethane in the strato- 

53 sphere (3x10^° s at 5 mbar) is significantly greater than one Jovian year (3.7x10^ s), 

54 whereas the acetylene lifetime (5x10'^ s at 5 mbar) is substantially shorter. Note that the 



55 solar cycle (3.4x10* s) which affects the photochemistry almost equals the Jovian year. 

56 In the presence of dynamical motions on seasonal timescales, the meridional abundance 

57 trend in a relatively short-lived species (such as acetylene) is expected to show significant 

58 change, whereas a long-lived species (ethane) should present a more constant distribu- 

59 tion. The logical means of testing this hypothesis is to intercompare the hydrocarbon 

60 abundances using the datasets of Voyager and Cassini, which encountered Jupiter at two 

61 different seasons: near northern fall equinox, and 1.75 jovian years later just after northern 

62 summer solstice. 

63 Such a comparison is the objective of this paper. We begin by analyzing the Voyager 

64 spectra, retrieving for the first time with a full radiative transfer modeling approach the 



65 meridional abundance trends of ethane and acetylene in 1979. We a. 



66 dataset, using a significantly revised spectroscopic atlas for ethane ( IVander Auwera et al. 



so reanalyze the CIRS 



67 120071 ) that was published since the first Cassini jovian hydrocarbons paper appeared (and 

68 is also used the IRIS data analysis here) . The intensities of this new list are revised upwards 

69 by some 44% relative to the GEISA 2003 list used in Paper I, resulting in a corresponding 

70 decrease in abundances. We then proceed to compare the abundance profiles at the two 

71 epochs (northern fall in 1979 versus summer in 2000) discussing the implications and 

72 drawing conclusions. 



2 Data Acquisition 



2.1 Cassini CIRS Observations 



75 The Casini CIRS instrument is a dual interferometer, with a mid-IR Michelson in- 

76 terferometer covering the range 6 00-1400 cm~ ^ , and a far-IR polarizing interferometer 



77 covering the range 10-600 cm ^ (JFlasar et al.l . 120041 ). The spectra analyzed here have 



78 been selected in the same way as in Paper I (described below), the only difference being 

79 that they are derived from a later version of the CIRS database. This has some im- 

80 provements to the calibration algorithm, including rejection of some 'bad' spectra (scan 

81 mechanism out of phase lock) and greater suppression of numerical artifacts from the 

82 Fourier transform. These mid-IR spectra were acquired at the highest spectral resolution 

83 of CIRS (0.48 cm~^ FWHM, Hamming apodized) during December 1-15 2000, and cover 

84 a latitude range of 70°S to 70°N with a spatial resolution of 6.0-3.5° great circle arc. The 

85 spectra were binned in 5° bins, and with auroral areas excluded as before. The effective 

86 emission angle 6 for the i^'^ bin was computed from the mean of the airmass (plane parallel 

87 atmosphere) of the Nj spectra as follows: 



6i = arcsec 



-— y^sec 

^^] j=l 



(1) 



89 For further details see Table I. 

[TABLE I appears here] 

90 2.2 Voyager IRIS Observations 

91 The two Voyager spacecraft passed much closer to Jupiter than Cassini (5.7 x 10^ km 

92 versus 9.8 x 10® km). The IRIS spectrometer carried on board was also a Michelson type - 

93 like the mid- infrared of CIRS - but with a single bolometer detector instead of two 1x10 

94 arrays. The IRIS spectra selected here are from the Voyager 1 inbound north/south map, 

95 taken over a 10.5 hour period beginning 2 days, 13 hours before closest approach at a 

96 range between 2.9 x 10® and 3.6 x 10® km. Nevertheless, the much larger field of view 

97 (FOV) of the IRIS bolometer (4.36 mrad) versus the CIRS mid-IR detectors (0.29 mrad 



Nixon et al.l. 



2009bi ) results in slightly larger projected field of view on the disk of 10 to 



99 12 degrees of latitude at disk center. The observation was made by repeatedly stepping 



100 the field of view from north to south along Jupiter's central meridian, with the rotation of 

101 Jupiter providing longitude coverage. The spectra were then averaged in 10° wide latitude 

102 bins, stepped every 5°. As the spectra were acquired at constant airmass for each latitude, 

103 the emission angles were simply averaged in each bin. See Table II for details. 

[TABLE II appears here] 

104 3 Data Analysis 

105 3.1 Model Atmosphere 



106 The model atmosphere used by us is identical to that of Pape r I. The initial tem pera 



19981 ). from 



107 ture profile was taken from the Galileo probe ASI measurements (iSeiff et al.l . 

lOB 4 bar to 0.4x10^^ bar. Hydro gen and helium have uniformly mixed vertic al abundances 



109 (0.863 and 0.134 respectively) ( iNiemann et al 



1998 



von Zahn et al. 



19981 ). Initial verti- 



110 cal abundance profiles for the three hydrocarbons (CII4, C2II6 and C2H2) were computed 

111 using a photochemical model, described in Paper I. Solar maximum conditions were used 

112 for the UV flux in this calculation for both the Cassini and Voyager epochs, as both 

113 encounters occured at similar points near the maximum in the solar cycle. 

114 The vertical profiles for the known constituents ammonia and phosphine were con- 
ns trolled by a three-parameter model: (1) deep abundance; (2) 'knee' pressure level, the 
lie transition level from deep, fixed, to varying abundance; (3) a fractional scale height used 
117 to decrease the abundance above the 'knee' level. The retrieval of ammonia and phos- 



118 phine parana e ters f r om the mid- in 



119 (jlrwin et al. 



200 



Abbas et al. 



rared spectrum has 



2004 



j een d i scussed extensively elsew here 



Fouchet et al. 



2004; 



Achterberg et al. 



20061). As 



120 these gases are not our focus here, we will remark only briefly on them concerning their 

121 effect on the hydrocarbon retrievals. 



With regard to PH3, the spectr al effect of the un b and (1121 cm ^) is significant in the 



region 900-1200 cm ^ (see Fig. 1 of 



Irwin et al. 



(I2OO4J )) and insignificant outside, including 



124 the hydrocarbon spectral ranges considered here. T herefore its parame ters were fixed at 



125 the values (6 x 10 ^, 1.0 bar, 0.3) recommended by 



Irwin et al 



mm . Regarding NH3, 

126 the spectral effect of the z/2 band at 950 cm~^ is most significant from ~800-1000 cm^^ 

127 (although the effects range from 750-1200 cm~^) , and so primarily affects the ethane 

128 spectral range. As it varies with latitude, it must be included in the retrieval. For our 

129 purposes the knee pressure was fixed at 0.7 bar and the both the deep pressure and 

130 fractional scale height were included as free parameters, with a priori values of (2.2 ± 

131 0.22) X 10^^ and 0.15 ± 0.05 respectively. Note also that ammonia is entirely separable 

132 from ethane due to the very different spectral signature. For discussion of the latitude 

133 and longitude variation of these retrieved parameters, see the above-referenced papers. 

134 Fig. [1] depicts the model atmosphere, which is divided into 71 vertical layers. 

Fig. 1. Appears Here. 

135 3.2 Forward Spectral Model and Retrieval Algorithm 

136 The computation of the s pectrum from the model atmosphere was carried out using 



137 the Nemesis computer code (llrwin et al 



20081 ). which also accomplishes the fitting and 



138 parameter retrieval. This code has been significantly validated in the past, having been 

139 successfully applied to model reflection and thermal infrared spectra of Venus, Mars, 

140 Jupiter, Saturn and Uranus; from a variety of spacecraft and ground-based facilities and 

141 instruments, including Galileo, Cassini, Mars Express, Venus Express, UIST. 

142 The Jovian model opacity is derived from two sources: (i) gas vibrational-rotational 

143 bands, and (ii) collisional-induced gas opacity. Jupiter's atmosphere is also known to 

144 contain stratospheric haze. The addition of haze (particulate) opacity was investigated 



8 



145 in paper I using the refractive index co-efficients of iKhare et al.l (J1984J ) or a simple grey 

146 absorber. However, the effect on the retrieval was not found to be statistically significant, 

147 i.e. the difference to the x^ was below the noise threshold. Therefore we have followed the 

148 approach of Paper I and do not include a haze opacity in our retrievals. The atlases for 

149 the gas bands are ident ical to Paper I, except for et hane, where we substituted the more 



150 recent, improvec 



151 ( Rothman et al. 



ist o 



Vander Auwera et al 



(120071 ) ■ now also available in HITRAN 200 8 



20091 ) and the 2009 update of GEISA (jJacquinet-Husson et al. 



2008 



). 



152 The collisio n-induced opacity uses the method of Borysow and co-workers (iBorysow et al. 



1985 . 



19881). 



154 The spectral calculation proceeds based on the correlated-A; approximation ( iGoody et al.l . 



19891), whereby tables of /^-coefficients are pretabulated in advance, at an appropriate 



range of pressure and temperature grid points covering the entire atmospheric range of 
interest. The spectral resolution was convolved into the /c-tables during the computation, 
so that separate tables were built for each gas at both the IRIS and CIRS resolutions: 
3.9 and 0.48 cm~ ^ FWH M of Hamming function instrument line shape (ILS) respectively 



160 see 



Nixon et al. 



2009al section 2.2 for further details of the apodization). This consti- 



tutes a change and improvement over the spectral modeling of Paper I, where the ILS 
was approximated by a triangle function, and is one source of the slight differences seen 
in the retrieved abundance of C2H2 from the Cassini CIRS spectra. See further discussion 
in ^EH 



165 Th e fitting of the spectrum to the data was achieved through the formalism of 



Rodgers 



166 (l2000l ). known as non-linear least squares optimal estimation. In essence, this is simi- 

167 lar to minimizing the reduced x^ difference between the model and data. However the 

168 goodness-of-fit measure, known as the cost function (p, additionally includes a constraint 

169 (or smoothing term) from a priori (assumed) information about the likely range of val- 

170 ues that the parameters can take. The cost function is therefore the quadratic sum of 



171 the spectral x^, plus a similar term that measures the deviation of the solution from the 

172 a priori co nstraint (see ^4.2 of Paper I). is minimized using the Levenburg- Mar quart 



173 technique ( jPress et al. 



19921 1 ■ along downhill gradients computed by taking the numerical 

174 partial derivatives of the solution vector (computed spectrum) with respect to the free 

175 parameters - also known as the functional derivatives. The retrieval terminates when a 

176 pre-set convergence limit is reached, in this case when the change in the cost function is 

177 less than 0.2% betwe en iterations. For further details regarding the Nemesis algorithm, 



178 see 



Irwin et al. 



fcoosh . 



179 3.3 Vertical sensitivity from functional derivatives 

180 Before commencing on the retrievals, we first used the forward model to investigate 

181 the vertical regions of information, by computing the functional derivatives 6I/6xij of 

182 the radiance / for each of the Xi parameters at each level j in the model atmosphere. The 

183 functional derivatives for the 0.48 cm~^ resolution Cassini CIRS spectra are described 

184 in Paper I, so we here describe the corresponding computations for the lower resolution 

185 (3.9 cm^^) Voyager IRIS spectra, noting differences from the CIRS sensitivities. 

Fig. 2. Appears Here. 

186 Figs. E] shows the functional derivatives for temperature in two spectral regions: in 

187 the hydrogen S(l) continuum at 600-650 cm~^ sensitive to the upper troposphere (~200 

188 mbar), and in the P and Q branches of the methane u^ band at 1260-1310 cm^^. The 

189 latter are sensitive to several altitude regions: 1-15 mbar in the middle stratosphere, but 

190 also at ~7 yubar in the upper stratosphere due to the optically thick center of the Q- 

191 branch at 1305 cm~^. Therefore, by fixing the abundances of H2 and CH4, we can sound 

192 the temperature structure at three different vertical levels. 

193 Figs. |3]-|4] show the functional derivatives for the gases C2H2 and C2H6. Acetylene has a 



10 



194 broad sensitivity of about 0.1-7 mbar in the stratosphere, but also has absorption at 300- 

195 800 mbar in the troposphere. Similarly, ethane has a range of 1-10 mbar the stratosphere, 

196 and 100-700 mbar in the troposphere. 

Fig. 3. Appears Here. 
Fig. 4. Appears Here. 

197 Fig. |5] summarizes the vertical sensitivity ranges of Voyager IRIS, and compares these 

198 values to those of Cassini CIRS, both calculated at the equator. See also Table III. The 

199 most important difference is that the hotband of C2H2 at ~731 cm~^, about 1 cm~^ in 

200 width, is resolved at the CIRS but not at the IRIS spectral resolution. Therefore, the 

201 elevated range of sensing for CIRS at 731 cm~^ (peak at 0.1 mbar, FWHM extending to 

202 0.01 mbar) relative to the z/5 Q-branch at 729 cm^^ (peak at 5 mbar, FWHM to 0.03 mbar) 

203 is lost. Fig. [6] and Table IV show similar information, but calculated at 58°N showing that 

204 the sensitivity ranges shift slightly with latitude. 

[TABLE III appears here] 

Fig. 5. Appears Here. 
[TABLE IV appears here] 

Fig. 6. Appears Here. 

205 Before leaving this topic, we also must consider whether the vertical retrieval range 

206 will be limited by the transition from local thermodynamic equilibrium (LTE) to non- 



207 LTE (NLTE) conditions. For the case of 



208 previously considered by 



Halthore et al. 



upite r's methane at least, the problem has been 



( 119941 ). who investigated the relaxation of dilute 



209 gaseous methane in bulk hydrogen. In that paper, the authors used a room temperature 

210 measurement for the collisional de-excitation time re = Ciq for the u^ band of methane 

211 (3.3 /im) in combination with a Landau- Teller formula to extrapolate the relaxation to the 

212 ground state under Jovian conditions. The same tq also applies to the z/4 band, as there 



11 



213 is fast transfer of vibrational energy (V-V) between these excitations. Their Equation 5 

214 may be re-arranged: 

215 re = 2.3175 x 10" V^ exp{AOT-^^^) (2) 

216 where p is pressure in Torr and T is temperature in K. In Fig. [7| we compare the ahitude- 

217 variation of re calculated for our initial equatorial atmosphere, versus th e spontaneous 



Halthore et al. 



(|l99J). We find that 



218 emission timescale Tp^ = A{q, taken as 0.4 s following 

219 rA ~ re at a pressure of ~ 0.8 fihai. We have also made an independent estimate of Aiq by 

220 averaging across the CH4 1/4 band in HITRAN from 1215-1315 cm~^, arriving at an even 

221 longer timescale of 1.7 s, putting the transition level well above the upper pressure limit 

222 our model atmosphere. The variation with latitude is small (using the profiles derived 

223 later in this paper). Therefore we conclude that our temperature retrievals are valid over 

224 the range indicated by the contribution functions, which probe pressures p > 0.8 /ibar. 

Fig. 7. Appears Here. 

225 4 Results 

226 The retrieval proceeded in the following two-step manner. Firstly, each latitudinal aver- 

227 age spectrum was analyzed to retrieve temperatures in the troposphere and stratosphere, 

228 using spectral portions of the hydrogen continuum (600-670 cm"^ and 760-800 cm"^) 

229 and methane 1^4 band (1225-1325 cm~^). Secondly, the temperature profile was fixed and 

230 the gas abundances of C2H2 and C2H6 allowed to vary, and their vertical distributions 

231 retrieved from their bands at 670-760 cm"^ (C2H2 u^) and 800-850 cm"^ (C2H6 I'g). 

232 Note that for reference, the retrieved 2-D temperature fields, gas mixing ratios for 

233 C2H2 and C2H6, and formal errors on these quantities, from both datasets (Voyager and 



234 Cassini epochs), have been archived on-line. See http://hdl.handle.net/2014/41479 



12 



235 and http://hdl.handle.net/2014/41480 For this reason, extensive tables of numeric 



236 quantities and errors are not enumerated in this paper, while some selected values are 

237 discussed in context. 

238 4-i Temperatures 

Fig. 8. Appears Here. 

239 Fig. [S] shows the retrieved temperatures for Voyager IRIS (top), Cassini CIRS (middle) 

240 and the difference (bottom). Considering first the stratosphere at 1 mbar, we see that the 

241 northern hemisphere has warmed considerably (~5 K) at this epoch, equivalent to 1.75 

242 Jovian years later. Overall, the Cassini lower stratosphere exhibits a more pronounced 

243 hemispheric asymmetry than the Voyager era, however the upper stratosphere (0.01 mbar) 

244 is colder (~ —10 K). Evidence for seasonal change is less apparent in the troposphere at 

245 100-400 mbar, as expected due to the much greater thermal inertia at these levels. 

246 Note th at the upper stratosphere region (0.01 mbar) was not considered in the earlier 



247 analysis of .Simon-Miller et al. 



( 1200 61 ). who also compare retrieved Jovian temperatures 

248 from Voyager and Cassini with an emphasis on derivation of wind fields. These authors 

249 used a different set of Cassini maps (from early January 2001) at lower (2.8 cm~^) spec- 

250 tral resolution than our selection (0.48 cm^^), although with higher spatial sampli ng in 



Simon-Miller et al. 



(12006f l. the 



251 latitude and longitude. Where our results overlap those of 

252 agreement is good (see their Fig. 2). 

253 We also see evidence for the the so-called jovian 'Quasi-Quadrennial Oscillation' (QQO), 

254 a cyclic temperature variation, most prominent in the lower stratosphere, that appears 

255 to observe an approximately four (terrestrial) year period, and may be analogous to the 

256 Earth's quasi-biennial oscillation (QBO, ~26 month period). See Fig. |H](c), at 20 mbar 

257 and 10°S, where a localized increase in temperature of ~8 K is apparent. The four-year 



13 



268 periodicity cannot be inferred from these two widely-spaced epochs alone, but when the 
259 time series of equ atorial stratospheric tempe ratures is plotted with higher temporal resolu- 



260 tion, e.g. Fig. 4 of 



Simon-Miller et al 



( 120061 ) , the pattern b ecomes apparent. The fou r- year 



261 cycle is not exact (hence 'quasi'): periodogram analysis^ 

262 cates two peaks with periodicities of 3.5 and 4.3 years. 



bz 



Simon-Mi 



Leovy et al 



leretaL 



fl2006[ ) indi- 



(Il99l[ ) have suggested 



263 that the QQO is caused by the stresses induced by vertically - propa g ating waves. For fur - 



264 ther details and discussion see 



Simon-Miller et al. 



(120061). 



Leovy et al. 



fll99lh : lFriedsonl(ll999[ l: 



Moses et al. 



(120041 ) 



266 4-^ Gas abundances 

267 The corresponding contour maps for the hydrocarbon abundances are shown in Figs. 

268 M and [ini Note that the vertical range of C2H2 derived from CIRS data has been cut to 

269 match the smaller vertical range of IRIS (see §3.3p . 

Fig. 9. Appears Here. 

270 The C2H2 difference plot (Fig. [9] (c)) shows the most dramatic difference, with the 

271 polar stratospheric abundances dropping at 1-10 mbar, and more intensely in the south 

272 than the north. This is also evident in Fig. |9] (b), where the upward-curved contours 

273 contrast markedly with the more level contours of Fig. |9] (a). 

Fig. 10. Appears Here. 

274 Ethane, in contrast, shows a much smaller change from Voyager to Cassini, although 

275 the Cassini era appears to have a more uniform stratospheric distribution, whereas the 

276 Voyager period is characterized by a positive bias towards the north. Also, the CIRS epoch 

277 shows systematically higher abundances over much of the atmosphere. 
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278 5 Discussion 

279 5. 1 Comparison to Paper I 

Fig. 11. Appears Here. 

280 The first reference point for this study is the previously published retrievals of hydro- 

281 carbon abundances from CIRS data (Paper I). The deviation between the present work 

282 and Paper I is plotted in Fig. [TTl For C2H2, the results are almost identical in regions 

283 where the information content is maximum, e.g. at 2 mbar in the equatorial stratosphere 

284 the change is just 1%, and at 300 mbar in the equatorial troposphere the change is 10%. 

285 The maximum change was at 62°S, 300 mbar where the deviation reached 25%. These 

286 differences are well within the overall error bars of ~40-50%, and may be attributed to 

287 small changes in the modeling method, such as the improved instrumental line shape (see 

288 §3.2p . See on-line data for full data and errors (§!]). 

289 For C2H6, the change at 5 mbar in the equatorial stratosphere was 40%, with a 20% 

290 change at 300 mbar. However, these differences are relative to smaller error bars of 15- 

291 20%, and are therefore clearly significant. These changes can be explained by the new line 

292 list used for C2H6 (see §3.2p . which has higher line intensities than used in Paper I, by 

293 some 44% over the whole band. 

294 5.2 Comparison to other previous studies 



295 In the recent review on the Jovian stratosphere by iMoses et al.l ( 120041 . §7.2.3 ), de- 

296 scribes the history of abundance measurements of C2H2 and C2H6 in the Jovian non- 
297 auroral stratosphere, spanning the period from 1973 to 1997. Their Table 7.1 lists all 
298 major observational measurements, totaling 19 published studies in all, broken down into 
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299 16 of C2H2 and 23 of C2H6. These varied researches cover Jupiter's stratosphere, mainly 

300 at low latitudes, but encompassing every season, and at pressures from 80 to 0.005 mbar. 

301 In light of the present findings, it is now clear that Jupiter's stratosphere is a dynamic 

302 environment that exhibits seasonal change in the abundances of the trace hydrocarbons 

303 at all latitudes and pressure levels. Therefore, direct comparison between our findings and 

304 these previous studies is not universally appropriate. Clearly, the most meaningful way to 
306 proceed is to compare results obtained at similar seasons to the results we have obtained 
306 for the Cas sini and Voyager ep ochs. Table V compares a relevant subset of observations. 



307 taken from 



Moses et al. 



(120041 . Table 7.1), that apply to the same season as either the 

308 Voyager or Cassini encounters. We have added error bars where those were available in 

309 the original sources. 

[TABLE V appears here] 



310 Relative to our Cassini results, we found 



311 only, which is the study of 



Livengood et al. 



only one applicable comparison, for C2H6 



( 119931 ) based on data from 12/89, about 1 



312 full Jovian year before Cassini. The Cassini-derived value of (9.4 ± 1.8) ppm at 1 mbar is 

313 somewhat greater their value of (3.6 ± 1.3) ppm, at 0.3-3 mbar, similar to our conclusion 

314 in Paper I. We note however that we would achieve agreement at their lower bound of 

315 3 mbar where we find an ethane VMR of (4.8 ±0.8) ppm. The other studies for C2H6 in the 

316 4-5 mbar region also agree with our results when the respective error bars are considered. 

317 Note that the entire set of past measurements for C2He is generally self-consistent, in 

318 marked contrast to those of C2H2, discussed hereafter. 



319 With re gard to C^H?, our results are compatible with those of 



320 mbar) and 



Noll et al. 



Owen et al 



(Il980h (10 



19861) ( 3 mb ar). We find about 1 orde r of rn agnitude less C2H2 



321 at 10 mbar thaniClarke et al. 



322 magnitude greater than 



(Il982h and 



Bezard et al. 



Gladstone and Yung) (119831 ). though 1 order of 



(jl995[ ). Note that the prior results at 10 mbar are 
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323 mutually inconsistent. The reason for this is unclear, but the lower abundance compared 

324 to C^Hfi may be pa rtly responsible. At 1 mbar, our result appears to be compatible with 



Bezard et al 



(I1995I ). if some error is allowed on those results. 



326 5.3 Chemistry and Dynamics 

327 The meridional distribution of a gas species depends upon the relative chemical and 

328 transport lifetimes. In the absence of meridional transport, for a species that is short- 

329 lived photochemically (compared to seasonal time-scales) we expect that the observed 

330 distributions in latitude should, to first order, follow the incident seasonal insolation. This 

331 was noted in Paper I for C2H2, whose lifetime is short (3 x 10^ s at 5 mbar) compared to 

332 both a Jovian season (1 x 10^ s), and the inferred tim escale for meridional mixing from 



333 analysis of Comet Shoemaker-Levy 9 debris (6 x 10^ s. 



Lellouch et al 



2OO2I) 



334 For species with photochemical lifetimes approaching seasonal time-scales, a phase lag 

335 in response to the insolation is expected with a diminished peak to peak variation in the 

336 mixing ratio. And finally for long lived species (e.g. ethane, with lifetime 3 x 10^*^ s at 

337 5 mbar) only small perturbations about annual global averages are expected. Meridional 

338 transport will wash out any differences. Additionally the vertical transport required to 

339 close the meridional cells will cause changes in the mixing ratios at various pressures 

340 depending upon the vertical profile of the mixing ratio, e.g. vertical transport bringing 

341 down C2H2 rich air to a given pressure level. 

342 Figure [12] compares the photochemical lifetimes (halving and doubling) for C2H2 and 

343 C2H6 from the 1-D model calculations described in Paper I to the Jovian year, the solar 

344 cycle, and the eddy mixing time (H^/K, where H is the mixed atmospheric scale height 

345 and K is the eddy diffusion coefficient). The acetylene lifetimes are less than the seasonal, 

346 solar cycle, and eddy mixing times, while the opposite is the case for ethane. In the 
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347 model both the acetylene and ethane mixing ratios increase with decreasing pressure, so 

348 downward transport of air parcels will increase their relative abundances and conversely 

349 for rising parcels. However, ethane is much more uniformly mixed with height due to 

350 its chemical lifetimes being longer than the eddy transport time. Thus vertical transport 

351 will have more of an effect on the C2H2 mixing ratio than for ethane, as enriched air is 

352 downwardly advected. 

Fig. 12. Appears Here. 

353 For acetylene and equinoctial conditions, we expect a more symmetric distribution 

354 compared to solsticial. Our results seem to validate this hypothesis, as the Voyager dis- 

355 tribution of C2H2 appears much more symmetric than the Cassini distribution, which 

356 is biased towards the northern hemisphere - as expected for northern fall equinox and 

357 summer solstice conditions respectively. For ethane, the reverse is seen: the Cassini (sol- 

358 stice) epoch shows a more uniform latitude distribution than the Voyager (equinox) epoch, 

359 which may indicate a phase lag between dynamics and seasons. 

360 The higher abundances of C2H6 seen at most latitudes by CIRS may in part be due to 

361 the increased mean insolation at northern solstice relative to fall equinox: Jupiter was 5% 

362 closer to the Sun in December 2000 compared to March 1979, resulting in 10% higher solar 

363 flux (21% flux variation possible due to solar distance over an entire year). In addition. 



364 we note that the increase in C2H6 tow ards the po 



es at b oth epochs may be partially due 



200l|). 



365 to particle-induced auroral chemistry (JYelle et al. 

366 5.4 Comparison to Saturn 



367 Our flndings for Jupiter closely parallel recent studies of the minor hydrocarbons in 

368 Saturn's atmosphere, which have mostl y found a decr e ase in acetylene and constant or 

369 increasing ethane from equator to pole. 



Howett et al. 



(120071 ) ■ using Cassini CIRS data. 



370 found a decrease in the abundance of C2H2 at 2 nibar from 30-68°S by a factor 1.8, 

371 although there was an initial increase from 18-30°S. Ethane on the other hand increased 



372 from the equator to the south pole by 2.5. 



Hesman et al 



( 120091 ). combining results from 



373 Cassini CIRS and ground based spectroscopy with the IRTF and Celeste, found a decrease 

374 in C2H2 by a factor 2.7 from 5-75°S, although there was a rise again at 87°S. Ethane was 

375 constant from equatorial to southern mid-latitudes, where it began to rise in abundance, 

376 finally doubling at the south pole. 

377 Both these studies u sed nadir-geometry spectra, limiting the range of vertical sensitiv- 



378 ity. 



Guerlet et al. 



(120091 ) however have analyzed Cassini CIRS limb spectra to yield a more 

379 detailed vertical picture of the hydrocarbon meridional abundance variations. Acetylene 

380 is found to be sharply peaked at the equator at 1 mbar, with a gradual decrease through 

381 southern and northern mid-latitudes, while ethane is nearly flat from 50°N-50°S, where 

382 is rises towards the south pole. 

383 Overall, we see some intriguing similarities between the meridional variations of the 

384 secondary hydrocarbons on Jupiter and Saturn, although the details differ. The implica- 

385 tion is that researches into understanding the origins of these distributions must proceed in 

386 parallel, and advances in understanding the mechanisms of one body are likely to improve 

387 our appreciation of the other. Note that on Saturn however we must be mindful of the 

388 effect of the ring shadow, which is seasonally varying over a 29.5 (earth) year period. This 

389 has the effect of shielding the low latitudes, reducing both the photolysis of methane, and 

390 the higher hydrocarbons, which must be accounted for correctly in photochmeical models. 
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391 6 Summary and Conclusions 

392 In this paper we liave analyzed two infrared spectral datasets of Jupiter: (i) the Voyager 

393 IRIS dataset, pertaining to the northern fall equinox of 1979, and (ii) the Cassini CIRS 

394 dataset, taken near the northern summer solstice of 2000. We have modeled the spectra 

395 to retrieve, firstly atmospheric temperatures through the z/4 band of CH4 at ~ 1305cm~^, 

396 and secondly the abundances of C2H6 and C2H2 based on the derived temperatures. Our 

397 results extend from the upper troposphere to the upper stratosphere for temperature 

398 and C2H2 abundance, and cut-off slightly lower in the middle stratosphere for C2H6. The 

399 results comprise the first comprehensive 2-D comparison of the abundances of C2H2 and 

400 C2H6 in Jupiter's atmosphere at two different seasons. The results have been archived 

401 on-line for reference, as described in ^ 

Our new findin gs confirm, improve and extend the results of our previous paper 



403 (JNixon et al. 



20071 ). The distribution of C2H2 at the Cassini epoch is confirmed, whilst 



404 that of CpiKfi is corrected, largely du e to revised spectral line intensities pubhshed in the 



405 interim ( IVander Auwera et al. 



20071 ). We have also extended the Cassini distributions to 



406 a second Jovian season (Voyager), showing important similarities and differences. 

407 For ethane, the 2-D distribution is quite similar, although we see that the Cassini 

408 picture is more symmetric than that of Voyager, where C2H6 is somewhat depleted in 

409 the north relative to the south. Also, the Cassini abundances are mostly higher overall, 

410 perhaps due to smaller solar distance. For acetylene the difference is much more dramatic, 

411 showing a uniform meridional profile at equinox, but strong depletions towards both poles 

412 at solstice. Therefore the solstice distribution seems to match the annual-mean picture 

413 that would be predicted by a 1-D photochemical model. 

414 It is clear that Jupiter's stratosphere is indeed a complex environment, where chemistry 
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415 is coupled to dynamics in ways that are poorly understood at present. Much additional 

416 modeling work is required before these seasonal changes can be confidently interpreted, 

417 and comprehensive radiative-dynamical-chemical models in two or three dimensions are 

418 urgently needed. The rewards of successful modeling of this data are large, and include 

419 the significant advancement of our understanding not only of the Jovian atmosphere, but 

420 also that of Saturn and the other outer planets as well. 
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TABLE I 
Cassini CIRS Observational Data 



Mean No. of FP3 No. of FP4 Mean Emis. NESR 

Latitude Spectra Spectra Angle (°) 760 cm~ 



NESR 
^ 1300 cm^i 



-66.30 


19 


-62.24 


73 


-57.28 


172 


-52.49 


231 


-47.54 


306 


-42.31 


559 


-37.36 


855 


-32.48 


1171 


-27.84 


644 


-22.32 


446 


-17.37 


427 


-12.20 


471 


-7.93 


437 


-2.65 


442 


2.50 


505 


7.55 


413 


12.56 


580 


17.52 


604 


22.44 


525 


27.33 


478 


32.47 


464 


37.51 


491 


42.56 


430 


47.37 


437 


52.29 


381 


57.16 


334 


62.24 


220 


67.14 


96 



7 

67 
139 
234 
313 
332 
366 
402 
377 
523 
501 
2144 
1149 
503 
551 
536 
558 
568 
492 
568 
548 
558 
489 
400 
402 
294 
187 

97 



71.21 
67.76 
64.13 
60.81 
56.93 
52.92 
48.96 
45.76 
38.73 
36.29 
32.23 
34.77 
27.51 
23.41 
21.89 
22.49 
26.31 
28.16 
29.66 
34.89 
38.67 
41.09 
46.05 
50.57 
55.60 
59.12 
63.54 
66.42 



3.50 X 10 
1.92 X 10 
1.50 X 10 
1.36 X 10 
1.09 X 10 
0.94 X 10 
0.94 X 10 
1.14 X 10 
1.07 X 10 



1.07 
1.02 
1.01 



10 
10 
10 



1.07 X 10 

1.07 X 10 
1.01 X 10 

1.08 X 10 

1.07 X 10 
1.03 X 10 

1.09 X 10 

1.08 X 10 

1.09 X 10 
1.17 X 10 

1.14 X 10 
1.16 X 10 

1.15 X 10 
1.23 X 10 
1.43 X 10 
1.94 X 10 



9 

9 

9 

9 

9 

9 

9 

-9 

-9 

9 

9 

9 

9 

-9 

-9 

-9 

9 

9 

9 

9 

9 

-9 

-9 

9 

9 

9 

9 

9 



8.50 X 10-1° 
3.30 X 10-1° 
2.30 X 10-1° 



1.89 
1.71 
1.71 
1.74 



10-1° 
10-1° 
10-1° 
10-1° 



1.71 X 10-1° 
1.71 X 10-1° 
1.93 X 10-1° 
1.65 X 10-1° 
1.16 X 10-1° 
1.59 X 10-1° 
2.07 X 10-1° 
1.83 X 10-1° 
1.88 X 10-1° 
2.02 X 10-1° 
2.05 X 10-1° 
2.20 X 10-1° 
2.52 X 10-1° 
2.39 X 10-1° 
2.11 X 10-1° 
2.13 X 10-1° 
2.22 X 10-1° 
2.30 X 10-1° 
2.51 X 10-1° 
3.10 X 10-1° 
3.79 X 10-1° 



N. B. this table, included here for compl etenes s, is identical to TABLE I of INixon et all 

(|2003). 



TABLE II 
Voyager IRIS Observational Data 



Mean Latitude No. of spectra Mean Emission Angle 



-59.97 


27 


-55.02 


32 


-49.49 


43 


-45.04 


55 


-40.73 


46 


-34.29 


52 


-30.34 


62 


-25.01 


53 


-19.49 


63 


-15.00 


68 


-9.59 


74 


-4.80 


88 


-0.34 


79 


4.79 


73 


9.72 


76 


14.82 


65 


20.02 


55 


25.19 


58 


30.01 


60 


35.05 


58 


39.91 


59 


44.04 


48 


50.05 


40 


54.35 


37 


58.63 


22 



mission An 


gle NESRt 


60.39 


9.62 


X 10-1° 


54.65 


8.84 


X 10-1° 


49.76 


7.62 


X 10-1° 


45.34 


6.74 


X 10-1° 


40.43 


7.37 


X 10-1° 


34.06 


6.93 


X 10-1° 


30.49 


6.35 


X 10-1° 


25.33 


6.87 


X 10-1° 


20.91 


6.30 


X 10-1° 


19.52 


6.06 


X 10-1° 


17.16 


5.81 


X 10-1° 


15.22 


5.33 


X 10-1° 


12.62 


5.63 


X 10-1° 


10.32 


5.85 


X 10-1° 


10.56 


5.74 


X 10-1° 


12.01 


6.20 


X 10-1° 


16.67 


6.74 


X 10-1° 


21.55 


6.57 


X 10-1° 


26.04 


6.45 


X 10-1° 


30.50 


6.57 


X 10-1° 


35.87 


6.51 


X 10-1° 


40.10 


7.22 


X 10-1° 


45.35 


7.91 


X 10-1° 


49.86 


8.22 


X 10-1° 


54.34 


1.07 


X 10-°9 



'^ NESR = Noise Equivalent Spectral Radiance (l-o") in units of W cm ^ sr i/cm i. 



TABLE III 

Comparison of Altitude Sensitivity: Equator 











Lower 


Upper 


Wavenumber 


Primary 


Retrieved 


Troposphere 


Stratosphere 


Stratosphere 


cm~^ 


Opacity 


Parameter 


Range (mbar) 


Range (mbar) 


Range (mbar) 



625 


H2 


Temp 


98-339 (206) 
106 339 (206) 


722 


C2H2 


VMR 


74-339 (178) 
266 748 (516) 


729 


C2H2 


VMR 


84-374 (219) 
262 772 (525) 


731 


C2H2 


VMR 


87-392 (242) 
262 772 (533) 


822 


C2H6 


VMR 


137-597 (356) 
73 679 (432) 


1250 


CH4 


Temp 




1304 


CH4 


Temp 




1306 


CH4 


Temp 





0.93-31.27 (9.80) 
0.12 7.21 ( 3.79) 

0.03-28.39 (4.59) 
0.17 7.45 ( 3.79) 

0.01- 2.83 (0.15) 
0.17 7.45 ( 3.85) 

0.70-13.97 (4.38) 
0.86 10.45 ( 4.38) 

1.59-21.59 (5.95) 
2.37-15.39 ( 6.14) 

0.92- 7.21 ( 3.27) 
1.01- 7.69 ( 3.61) 

0.40- 2.79 (1.25) 
0.76 9.96 ( 3.12) 



0.003-0.025 (0.006) 
0.003 0.019 (0.007) 



VMR = Volume Mixing Ratio. 

Italic text is Cassini CIRS; bold text is Voyager IRIS. 

In each case the range and peak (in parentheses) of the sensitivity is given. 



TABLE IV 

Comparison of Altitude Sensitivity: North 











Lower 


Upper 


Wavenumber 


Primary 


Retrieved 


Troposphere 


Stratosphere 


Stratosphere 


cm~^ 


Opacity 


Parameter 


Range (mbar) 


Range (mbar) 


Range (mbar) 



625 


H2 


Temp 


63-254 (159) 
66 262 (156) 


722 


C2H2 


VMR 


71-284 (106) 
206-587 (399) 


729 


C2H2 


VMR 


74-323 (111) 
202-597 (399) 


731 


C2H2 


VMR 


73-392 (242) 
202-597 (399) 


822 


C2H6 


VMR 


90-525 (318) 
135 551 (345) 


1250 


CH4 


Temp 




1304 


CH4 


Temp 




1306 


CH4 


Temp 





3.22-39.19 (16.95) 
1.51 11.15 ( 5.31) 

0.09-39.19 (17.23) 
1.54 11.89 ( 5.48) 

0.06- 2.83 (0.16) 
1.54 12.08 ( 5.57) 

1.09-19.91 (2.83) 
2.70 21.59 ( 9.04) 

2.09-25.77 (9.64) 
2.37 21.59 ( 7.95) 

1.19- 6.44 (2.97) 
0.77- 6.99 ( 2.88) 

0.39- 2.45 (1.31) 
0.59 14.43 ( 2.30) 



0.002-0.007 (0.004) 
0.003 0.015 (0.005) 



VMR = Volume Mixing Ratio. 

Italic text is Cassini CIRS; bold text is Voyager IRIS. 

In each case the range and peak (in parentheses) of the sensitivity is given. 



TABLE V 

Comparison to Previous Results 



Gas Ref. Press. Lat.'^ Date Mole Voyager 79 Cassini 2000 

Species (mbar) Frac. Mole Frac. Mole Frac. 



C2H6 


FBI 


0.005 


16.5°N 


07/79 


(2.5 ±1.7) ppm 


N/A — 




G83 


4 


30°S-30°N 


12/78 
& 6/79 


(6.6 ± 5.3) ppm 


(2.3 ± 0.7) ppm — 




N86 


5 


8°S-8°N 


06/80 


(2.6 ± 0.9) ppm 


(1.9 ±0.6) ppm — 




L93 


0.3-3 


45-55° N 


12/89 


(3.6 ± 0.9) ppm 


— (9.4±1.8) ppm* 


C2H2 


O80 


10 


30°S-30°N 


1/79 
&4/79 


19ppb 


(13 ± 5) ppb — 




C82 


10 


30°S-30°N 


1978-1980 


(100 ± 50) ppb 


(13 ± 5) ppb — 




G83 


10 


30°S-30°N 


12/78 
& 6/79 


(100 ± 10) ppb 


(13 ± 5) ppb — 




W85 


30 


30°S-30°N 


1978-1980 


(30 ± 10) ppb 


N/A — 




N86 


3 


8°S-8°N 


6/80 


(34 ± 17) ppb 


(62 ± 20) ppb — 




B95 


0.1 


10-36° S 


1/91 


4 ppm 


N/A — 




B95 


1 


10-36° S 


1/91 


100 ppb 


(250 ± 100) ppb — 




B95 


10 


10-36° S 


1/91 


1 ppb 


(16 ± 6) ppb — 



^Voyager/Cassini data calc ulated at mean lati tude of range, where a pplica ble. -I-At 3 

mbar. Re f erence s: O80 =IOwenet all (flMI: F81 = IFestou et al.1 (flM); C82 = 

IClarke et all (119821 ) ; G83 = IGlads tone and Yund ffT983l); W85 = IWagener et aH ffTOSST ): 



N86 = INollet all (fl986[ ): L93 = .Livcngood'd^lH (|l99l l: B95 = IBezard et ^ ( fl995[ ) . 



572 Figure Captions 

573 Fig. [1] Initial assumed a priori temperature (left) and gas profiles (right) for Jupiter 

574 atmospheric model. 

575 Fig. [2] Functional derivatives for temperature from the Voyager IRIS spectral data: (a) 

576 from the H2 continuum region giving tropospheric information, and (b) from the CH4 z/4 

577 band, sensitive to the lower stratosphere near 5 mbar and the upper stratosphere at 5 

578 /ibar. 

579 Fig. E] Functional derivative for gas abundance from the C2H2 1^5 band. Note the 

580 emission emanating from the stratosphere at 1 mbar, and also tropospheric absorption at 

581 ~200 mbar. Some information is available from the upper stratosphere at 10-100 yubar 

582 due to the Q-branch (730 cm"^). 

583 Fig. m Functional derivative for gas abundance from the C2H6 i^g band. Note the 

584 emission emanating from the stratosphere at 5 mbar, and also tropospheric absorption at 

585 ~200 mbar. 

586 Fig. [5] Examples of the FWHM sensitivity (symbols at bar ends) and peak sensitivity 

587 for the contribution functions at various representative wavenumbers, comparing Cassini 

588 (x) to Voyager (o), computed for Jupiter's equator. Even though the spectral resolutions 

589 are quite different (0.5 cm^^versus 4.3 cm^^), the vertical information regions are similar 

590 except for the C2H2 hotband near 731 cm~^ for which the peak has moved considerably 

591 lower for IRIS due to lower spectral resolution. 

592 Fig. [6] As in previous figure, but computed for 58°N. 

593 Fig. [7] Comparison of the vertically- varying coUisional relaxation time (t c.) for the 1^4 



594 band of CH4 to the spontaneous emission timescale ta, assumed 0.4 s following 



Halthore et al. 



595 (I1994I ). The break-down of LTE conditions is computed to be at ~ 0.8 fihai. 
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596 Fig. [8] Retrieved temperatures from Voyager IRIS data (a) and Cassini CIRS (b). (c) 

597 shows the difference: Cassini- Voyager. This may be compared to Simon-Miller et al. (2006) 

598 Fig. 2. We find that the stratosphere at ~2 mbar is warmer in the north in 2000 compared 

599 to 1979, and at 0.01 mbar is cooler. 

600 Fig. [9] Comparison of retrieved abundances for C2II2 at the (a) Voyager and (b) Cassini 

601 epochs. The vertical range has been restricted to match the smaller span (IRIS) (see Figs. 

602 Oandin]). Abundances in the lower stratosphere (~ 10 mbar) have dropped substantially 

603 in 2000 at the poles relative to 1979. 

604 Fig. (TD] Comparison of retrieved abundances for C2H6 at the (a) Voyager and (b) 

605 Cassini epochs. The changes seen are much less dramatic than for C2H2, nevertheless 

606 the Cassini (2000) stratospheric abundances appear to be higher and more symmetric in 

607 latitude compared to the Voyager (1979) contours, which show a positive bias to southern 

608 hemisphere. 

609 Fig. [TT] Comparison of the retrieved abundances of C2II2 and C2H6 from the Cassini 

610 CIRS dataset, taking the difference between the 2007 and 2010 retrievals. The differences 

611 for acetylene, due mainly to changes in the instrumental lineshape assumed (see text 

612 for discussion), are typically 1-10% across most of the range (c.f. formal retrieval errors 

613 of ~ 40 50%), reaching maximum divergences of +20% and -40% in the northern 

614 stratosphere and troposphere respectively, which approach the size of the inferred random 

615 errors. For ethane, the differences are greater, with large changes (-40 to -60%) occuring 

616 in the mid-stratosphere at 10 mbar at the peak of the contribution function, and are 

617 significantly larger than the random errors (15-20%). This is due to the revised z/g band 

618 intensities used for C2II6, as discussed in §3.2. 

619 Fig. [T2] Photochemical production (solid color lines) and loss (broken color lines) times 

620 (time for mixing ratio of species to increased/decreased by a factor of two) for C2H2 (red) 
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621 and C2H6 (blue) for the model used in iNixon et al.l (120071 ) . These are compared to the 

622 eddy mixing time, the Jovian year, and half the period of the solar cycle. 
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(a) CjH^: Voyager/IRIS VMR 
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a) CjH,: Voyager/IRIS VMR 
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(c) CjH^: CIRS-IRIS difference in log WR 
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H, VK/R: CIRS 2010-2007 Re:rieval Difference % 
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